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FOREWORD 

This  report  describes  a  portion  of  the  results  obtained  on  NASA  Grant  NSG 
3044.  This  work  was  done  under  subcontract  to  the  University  of  Illinois, 
Urbana,  with  Prof.  S.S.  Wang  as  the  Principal  Investigator.  The  prime  gran¬ 
tee  was  the  Massachusetts  Institute  of  Technology,  with  Prof.  F.J.  McGarry  as 
the  Principal  Investigator  and  Dr.  J.F.  Mandell  as  a  major  participant.  The 
NASA  -  LeRC  Project  Manager  was  Dr.  C.C.  Chamis . 

Efforts  in  this  project  are  primarily  directed  towards  the  development  for 
finite  element  analyses  for  the  study  of  flaw  growth  and  fracture  of  fiber 
composites.  This  report  presents  a  convenient  method  of  analysis  for  cracks 
emanating  from  holes  in  anisotropic  materials. 
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ABSTRACT 

A  method  of  analysis  has  been  developed  for  cracks  emanating  from  a  cir¬ 
cular  hole  in  an  off-axis  unidirectional  fiber-reinforced  composite.  This 
analysis  method  is  formulated  on  the  basis  of  conservation  laws  of  elasticity 
and  of  fundamental  relationships  in  anisotropic  fracture  mechanics.  The  prob¬ 
lem  is  eventually  reduced  to  a  system  of  linear  algebraic  equations  in  mixed¬ 
mode  stress  intensity  factors.  Superiority  of  this  current  method  to  other 
approaches  in  investigating  the  problem  with  a  very  complicated  crack  geometry 
and  material  anisotropy  is  demonstrated  when  used  in  conjunction  with  any 
numerical  method  such  as  a  finite  element  analysis.  This  method  is  found 
to  be  convenient  and  accurate.  Mixed-mode  stress  intensity  factors  and  the 
associated  energy  release  rates  in  the  crack  problem  are  determined  for  the 
composites  with  various  fiber  orientations.  Solutions  for  both  single  and 
double  cracks  emanating  from  the  edge  of  a  hole  in  the  composites  are  pre¬ 
sented  also  to  illustrate  the  fundamental  nature  of  the  problem. 


1.  INTRODUCTION 


Fracture  of  fiber-reinforced  composite  materials  containing  stress  con¬ 
centrations  such  as  cutouts  has  been  of  great  concern  for  years  in  the  design 
and  analysis  of  composite  structures.  Composites  with  circular  cutouts  of 
various  sizes  exhibit  the  so-called  hole-size  effects,  which  have  recently 
received  considerable  attention  from  researchers  because  of  their  significant 
variations  of  material  strength  reduction  [1],  Thus,  the  problem  of  crack  exten¬ 
sion  from  a  hole  in  composites  has  been  subjected  to  extensive  investigation 
both  analytically  and  experimentally.  While  the  importance  of  the  problem  has 
been  well  recognized,  understanding  of  the  fundamental  nature  of  the  problem 
is  still  very  limited  due  to  its  inherent  complexities.  Various  studies  have 
been  conducted  to  examine  the  stress  concentration  around  cutouts  in  infinite 
and  finite-dimensional  anisotropic  composites  and  have  achieved  significant 
success  [2-5],  However,  research  progress  in  studying  crack  initiation  and 
growth  from  a  hole  in  anisotropic  composite  has  been  relatively  slow.  The 
difficulties  involve  not  only  the  complex  crack  geometry  but  also  the  strong 
material  anisotropy.  In  addition,  the  crack-tip  singular  behavior  and  its 
interaction  with  the  hole  and  the  finite-dimensional  boundaries  of  the  com¬ 
posites  are  inherently  associated  with  the  problem.  These  complications  make 

it  extremely  difficult  to  be  solved  mathematically,  even  for  a  purely  elastic 
case , 

Early  work  by  Bowie  [6]  has  provided  solutions  for  the  problem  of  a  crack 
(or  cracks)  emanating  from  a  circular  hole  in  an  isotropic  infinite  plate  sub¬ 
jected  to  uniaxial  and  biaxial  stresses.  The  stress  intensity  factor  coeffi¬ 
cients,  expressible  as  a  function  of  the  ratio  of  crack  length  to  hole  radius, 
are  given  explicitly  in  [7]  in  a  tabulated  form.  Problems  of  crack  extension 
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from  a  hole  in  a  finite-dimensional  isotropic  strip  have  also  been  investi¬ 
gated  by  many  researchers  using  various  numerical  and  experimental  techniques 
[8-10],  However,  relatively  few  studies  of  the  problem  for  anisotropic  com¬ 
posites  have  been  reported.  Waddoups  et  al  [1]  were  among  the  first  to 
examine  the  failure  of  finite-dimensional  composites  with  circular  cutouts 
using  a  fracture  mechanics  approach.  Although  Bowie's  solution  for  an  iso¬ 
tropic  infinite  solid  was  assumed  to  be  held  for  the  finite-dimensional 
anisotropic  composites  in  the  study,  results  seem  to  predict  the  proper  trend 
of  the  fracture  behavior  of  the  composites.  Recently,  a  number  of  alternative 
approaches  have  been  proposed  [11-13].  Because  of  the  complex  nature  of  the 
problem  and  the  aforementioned  difficulties  involved,  most  of  the  solutions  are 
approximate  in  nature.  While  some  of  the  approximate  studies  provide  reasonable 
correlation  between  solutions  and  experimental  data,  no  vigorous  analysis,  to 
the  authors'  knowledge,  of  the  current  problem  has  been  reported. 

This  paper  describes  an  analytical  method  for  crack  extension  from 
a  hole  in  a  unidirectional  fiber-reinforced  composite  strip.  The  unidirectional 
fiber  composite,  which  has  the  simplest  laminate  configuration,  is  deliberately 
chosen  in  this  study  to  provide  basic  understanding  of  the  problem  and  to  serve 
as  a  reference  for  further  investigation  of  more  complex  multilayered  composites. 
General  concepts  underlying  the  basic  formulation  of  the  current  approach  are 
discussed  in  the  next  section.  Solution  procedure  for  numerical  computation 
is  presented  in  section  3.  The  mixed-mode  stress  intensity  factors  character¬ 
izing  the  crack-tip  behavior  are  determined  explicitly  in  section  4.  Effects 
of  fiber  orientation  on  deformation  and  crack  extension  are  examined  for  both 
single  and  double  crack  cases.  Of  particular  interest  are  the  energy  release  rates 
during  crack  extension;  their  values'  are  evaluated  for  several  cases  to  illustrate 


the  fundamental  nature  of  the  crack  problem.  Crack  initiation  and  growth 
from  a  circular  cutout  in  a  multilayered  fiber  composite  give  an  extremely 
complicated  three-dimensional  crack  geometry  and  have  to  be  studied  by  a 
fully  three-dimensional  fracture  analysis.  Some  progress  has  been  made  by 
the  authors  and  will  be  reported  elsewhere  [14]. 
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2.  FORMULATION 

The  current  approach,  using  conservation  laws  of  anisotropic  elasticity 
and  fundamental  relationships  of  anisotropic  fracture  mechanics,  can  provide 
rigorous  solutions  for  the  mixed-mode  fracture  parameters  without  the  need  of 
solving  the  corresponding  complex  boundary  value  problem.  The  method  has  been 
demonstrated  to  be  superior  to  other  approaches  in  its  mathematical  rigorous¬ 
ness  and  operational  simplicity  [15],  since  the  conservation  integral  can  be 
evaluated  conveniently  along  a  path  removed  from  the  crack-tip  region.  Details 
of  the  formulation  for  the  mixed-mode  crack  problem  in  anisotropic  composites 
have  been  reported  elsewhere  [16];  only  a  brief  description  of  the  fundamentals 
of  the  analysis  is  given  in  this  section. 

Consider  a  unidirectional  composite  plate  with  rectilinearly  aniso¬ 
tropic  elastic  properties.  Deformation  and  stress  fields  in  the  solid  are 
governed  by  the  generalized  Hooke's  law, 

£i  =  ^aijaj,  ij  =  1,2,6  (1) 

j 

where  a.,  are  compliance  coefficients  with  a..  =  a...  Compatibility  and 
ij  ij  Ji 

equilibrium  conditions  of  the  problem  lead  to  a  fourth  order  governing  partial 
differential  equation  with  a  characteristic  equation  of  the  form 

all  y4  ‘  2a16  y3  +  (2a12  +  a66)y2  "  2a26  p  +  a22  =  0  (2) 

Lekhnitskii  showed  [17]  that  the  roots  of  Eq.  2  are  always  complex  or  purely 
imaginary  and  occur  in  conjugate  pairs  as  p^ ,  p^  and  p^,  ^or  inf initesimal 

deformations  of  anisotropic  elastic  bodies,  a  form  of  general  conservation  laws 
for  the  plane  elasticity  problem  may  be  written  as 
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J. 

1 


3u. 


W  n.  -  a  n 

i  ik  k  3x. 

J  i 


ds  =  0 


(3) 


in  which  S  is  an  arbitrary  closed  boundary,  enclosing  a  portion  of  the  continuum 
with  unit  outward  normals  n_^,  and  W  is  the  strain  energy  density.  This 
conservation  law  can  be  reduced  [l8l  to  the  well-known  path-independent  J- 
integral  for  an  arbitrary  path  T,  which  begins  on  one  crack  surface,  encloses 
the  crack  tip,  and  terminates  on  the  opposite  face,  as  shown  in  Fig.  1, 


J  =  Jl  = 


9u. 


Wdy  -  T. 


1 


i  3x 


ds 


(4a) 


where 


T.  =  a . .  n. 
i  ij  J 


(4b) 


The  crack-tip  stresses  have  been  proved  to  possess  a  classical  fracture 

mechanics  r~i  singularity  in  a  linear  elastic  anisotropic  body  [19 ]  and  can 

be  characterized  by  stress  intensity  factors  K^.  and  (opening  and  in-plane 

shearing  modes,  respectively).  For  a  mixed— mode  fracture  problem,  the  correspon¬ 
ding  energy  release  rates,  Gj  and  G-^.,  may  be  related  to  K-j.  and  KI];  by 


GI  =  - 


K*  T 

T  a22  Im 


Vyi  +  V  +  Kn 


yl  y2 


(5) 


K 


JII 


~  a11  ImjK].I(p1  +  y2)  +  KI  ^  p2 


l 


(6) 


and 


Total 


GI  +  GII 


(7) 
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Using  Eqs .  5-7,  J  can  be  shown  to  be  related  to  the  stress  intensity  factors 
of  a  cracked,  anisotropic  composite  plate  subjected  to  in-plane  loading  by 

J  =  “ll  KI  +  al2  KI  KII  +  a22  KII  (8) 


where 


*22  -y  yl  +  y2 


a22  =  ~2~  Im^i  + 


yl  y2 


+  ~Y' 


It  should  be  noted  that  Eq .  8  alone  does  not  provide  adequate  information  to 
determine  and  separately  for  the  present  mixed-mode  crack  problem. 

If,  however,  the  conservation  law  for  two  independent  equilibrium  states 
[18]  is  introduced,  additional  information  may  be  obtained  to  remove  this 
deficiency.  Denoting  the  two  equilibrium  states,  1°  and  2°,  by  the  super¬ 
scripts  of  field  variables,  the  J-integral  for  the  superimposed  state  0°  may 


be  expressed  by 


j(0)  = 


where 


,(1,2) 


W<X’2)  dy 


m  3uf2)  .  SuP 

T  f 1  — ~ —  +  lf2)  — ^ —  ds 

x  3x  l  9x 
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r(k) 


w(k)  dy 


.  3ufk) 

I,(k)  -f-  is 

1  dX 


k  =  0,1,2  (14) 

(no  summation  on  k) 


„<»>  -  +  uf2) 

111 


(1  2) 

and  W  *  is  the  mutual  potential  energy  density,  denoted  by 


(15) 


(1,2)  (1)  (2)  (2)  (1) 

ijk£  i,j  k£ij  i,j  ^,5, 


(16) 


Recalling  the  J-K  relationship  in  Eq.  8  and  the  orthogonality  of  the 

(1  2) 

stress  intensity  factors,  M  ’  can  be  shown  [16]  to  have  the  form 


M(1,2)  =  2an  K^1}  kJ2)  +  K^2)  +  K^2)  K^} 


,  „(1)  f(2) 

2  22  KII  KII  (17) 


The  M-integral  in  Eqs .  13  and  17  deals  with  interaction  terms  only  and  can  be 
used  directly  to  determine  the  mixed-mode  stress  intensity  solutions  of  the 
present  problem. 


8 


3.  SOLUTION  PROCEDURE 


Equation  13  together  with  Eq.  17  provides,  in  fact,  sufficient  information 
for  determining  and  for  the  current  mixed-mode  anisotropic  fracture 

problem  if  proper  auxiliary  solutions  are  provided.  Introducing  the  first 
auxiliary  solution,  denoted  by  the  superscript  2a,  for  the  cracked  body  sub¬ 
jected  to  mode  I  deformation  only  so  that 


K 


(2a) 


=  1 


and 


K^a)  =  0 


(18) 


Eq.17  becomes 


(1,2a)  2a  +  a 

11  I  12  II 


(19) 


where  the  integral  M^’2a^  has  the  form 


M 


(1,2a)  . 


„  (1)  (2a)  , 

ijk£  Ui,j  Uk,X  y 


m  8u?2a)  .  .  9u^1) 

T (1)  i  .  T (2a)  l 

i  9x  i  9x 


ds 


(20) 


and  uf^  in  Eq.  20  can  be  determined  easily  along  a  properly  chosen 
integration  path  1.  Any  convenient  method  such  as  the  commonly  used  finite- 
element  analysis,  which  can  provide  accurate  far  field  solutions,  may  be  used. 

(  2  3.)  (2ci) 

u^  and  can  be  calculated  from  the  following  relationships: 


(2a) 

u .  = 

l 

A 

/  7T 

(21) 

and 

T (2a)  = 
i 

(2a) 

a.  .  n. 

ij  3 

(22) 

where 

a<2a)  =  - 

l 

/2irr 

gij^  ^i>1J2’e') 

(23) 

and  f^^  and  gf ^  are  given  explicitly  in  Ref.  [19], 
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Similarly,  a  second  auxiliary  solution,  denoted  by  the  superscript  2b,  for 
pure  mode  II  deformation  of  the  cracked  body  is  introduced  also  so  that 

K^2b)  =  0  and  K^2b)  =  1  (24) 


This  gives 


M^’2b)  =  a12  K™ 


2a  K(1) 
2  22  KII 


(25) 


where 


M 


(1,2b)  _ 


C.  uf1^ 
ijk£  i,j 


(2b)  , 
“k,i  dy 


a  (2b)  .  (1) 

TU)  i  ,  T (2b)  x 
i  3x  i  3x 


ds 


(26) 


The  quantities  T^2b^ 


/  4t  \ 

and  u^  in  Eq.  26  may  be  determined  by 

T<2b)  -  of2b> 

1  1J  J 


(27) 


where 


„<2b> 

ij 

l  di) 

(Ul,.i2 ,9) 

(28) 

and 

u(2b) 

1 

-  yf  4"V<v> 

(29) 

Thus,  Eqs .  19  and  25  provide  a  system  of  linear  algebraic  equations,  and 
and  can  be  determined  easily.  It  should  be  noted  that  in  the  process  of 

solving  for  and  ,  the  integrals,  M^b,2a^  and  M^b,2b^,  need  to  be 

evaluated  accurately  and  explicitly.  For  a  given  crack  geometry  and  loading 
condition,  this  can  be  achieved  by  integrating  the  functionals  in  Eqs.  20  and 
26  numerically  along  a  convenient  path  in  the  far  field  so  as  to  avoid  the 
crack-tip  complications .  In  this  paper  a  finite  element  analysis  is  used  in 
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conjunction  with  a  second-order  Gaussian  quadrature  for  evaluating  the  M-integral 
in  each  element  (Fig.  2)  by 


,d,2) 


N  3 

=  y  V  H  F  (x  ,y  )  As 
L.  ^  m  n  m  m  n 
n=l  m=l 


(30a) 


and  H  is  a  weight  coefficient  in  the  Gaussian  quadrature;  (x  ,y  ) ,  the 
m  mm 

Cartesian  coordinates  of  Gaussian  stations;  N,  the  total  number  of  elements 
associated  with  the  integration  path  F,  and  As^  is  a  segment  of  F  in  the  n-th 
element . 
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4.  RESULTS  AND  DISCUSSION 


The  study  of  mixed-mode  crack  extension  from  a  circular  hole  in  an 
anisotropic  fiber  composite  strip  has  been  conducted  successfully  using  the 


method  of  analysis  mentioned  in  the  previous  section.  The  traction  and  displace- 

fields 
d,2b) 


ment  fields,  T^  and  u^  ,  involved  in  Eqs.  20  and  26  for  determining 


and  M 


may  be  evaluated  along  a  properly  selected  path  removed  from  the 


crack  tip  by  a  finite-element  method,  which  employs  eight-node  isoparametric 
elements  of  the  displacement  model.  Accuracy  of  current  results  is  illustrated 
first  by  examining  a  crack— hole  problem  whose  solution  is  available  in  the 
open  literature.  Mixed-mode  crack-tip  stress  intensity  solutions  and  the 
associated  energy  release  rates  during  crack  extension  are  obtained  then 
for  various  cases  of  interest. 


4.1  Accuracy  of  Solutions 

Accuracy  of  current  solutions  for  the  mixed-mode  fracture  mechanics 
problem  has  been  studied  by  examining  the  problem  of  a  crack  emanating  from  a 
circular  hole  in  an  isotropic  finite— dimensional  plate  subjected  to  a  uniaxial 
in— plane  stress .  Figure  3  shows  the  geometry  and  individual  dimensions  of 
the  problem.  A  crack  of  length  a  is  assumed  to  emanate  radically  at  a  position  (R,<J>) 
on  the  hole  boundary,  where  R  is  the  radius  of  the  hole  and  <j>  is  the  direction 
along  which  the  crack  extends.  For  comparative  purposes,  the  problem  with 
<t>  =  0°  was  selected  for  evaluating  the  accuracy  of  the  solution  for  the  crack- 
hole  problem.  The  analysis,  employing  only  56  elements  in  the  numerical  com¬ 
putation,  gives  the  result  shown  in  Table  1,  which  is  in  excellent  agreement 
with  the  existing  solution  reported  by  Stalk  and  Orringer  [20],  A  detailed 
study  of  the  convergence  and  accuracy  of  the  present  method  of  analysis 
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for  mixed-mode  fracture  problems  has  been  reported  elsewhere  [15,16] .  The  results 
indicate  that,  in  general,  solution  accuracy  within  one  percent  deviation  from 
the  reference  ones  obtainable  in  literature  can  be  achieved  easily  using  only 
one-third  of  elements  in  other  highly  sophisticated,  singular  finite-element 
analyses  [20,21],  The  accuracy  of  the  current  solutions  is  further  reported 
to  be  insensitive  to  the  degree  of  discretization  in  the  numerical  calculation. 
This  situation  is  expected,  since  the  present  approach  bypasses  the  highly 
sensitive  and  complicated  crack-tip  region  and  evaluates  the  conservation 
integral  along  a  selected  path  in  the  far  field,  where  field  variables  are 
generally  smooth. 

4.2  Cracks  Emanating  from  a  Circular  Hole  in  Composites  with  Various  Fiber 

Orientations 

More  complex  fracture  problems  of  a  crack  (or  cracks)  emanating  from  a 
circular  hole  in  off-axis  unidirectional  fiber  composite  strips  are  examined 
in  this  section.  Since  the  response  of  the  problem  is  antisymmetric  in  nature, 
two  cracks  of  equal  length,  a,  are  assumed  to  be  initiated  from  the  boundary  of 
the  hole  in  the  composite , where  high  stress  concentration  occurs.  The  corre¬ 
sponding  single  crack  case  is  also  studied  for  comparative  purposes .  Figure  3 
shows  the  geometry  of  the  cracks  and  the  circular  hole  located  at  the  center  of 
a  finite-dimensional  2b  x  2L  anisotropic  panel  subjected  to  uniaxial  nominal 
stress,  a  ,  along  the  y-direction.  The  material  axes,  denoted  by  1  and  2,  of 
the  plate  are  oriented  with  respect  to  the  structural  axes  x  and  y  by  an  angle  0, 
as  indicated  in  the  figure.  The  local  stress  field  around  the  hole  (without  a 
crack)  in  the  finite-dimensional  composite  strip  is  solved  first  to  determine 
accurately  local  stress  concentrations  by  an  anisotropic  finite-element  analysis 
developed  in  the  present  study.  The  Tsai-Hill  stress  failure  criterion  [22 ]  is 
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introduced  then  into  the  analysis  to  determine  the  locations,  (R,<}>)  and  (R,Tr+<j>), 
of  crack  initiation  from  the  edge  of  the  hole  by 


F(ct. .) 
ij 


K] 

2 

°1  °2  A 

M 

UJ 

-  X2  + 

(Yj 

12  _ 
S2  “ 


(31) 


where  subscripts  1  and  2  refer  to  the  stresses  in  the  material  coordinates, 
and  X,  Y,  and  S  are  the  longitudinal,  transverse,  and  shear  strength  of  the 
unidirectional  composite,  respectively.  While  other  criteria  can  be  imple¬ 
mented  easily  in  the  current  analysis,  no  attempt  has  been  made  to  evaluate 
the  accuracy  of  different  failure  theories.  For  an  E-glass/epoxy  composite 
strip  panel,  X,  Y,  and  S  have  the  following  values  [23]: 

X  =  150  x  103  psi 

Y  =  4  x  10 3  psi 

and  S  =  6  x  10 3  psi 

Values  of  <j>  for  various  degrees  of  material  anisotropy  (i.e.  fiber  orientation) 
of  the  composite  are  obtained  by  Eq.  31  in  the  present  analysis.  Figure  4 
gives  the  location  of  crack  emanation  <j>  as  a  function  of  0  for  the  E-glass/ 
epoxy  composite  indicating  that  <j>  decreases  appreciably  with  increasing  0. 

It  is  a  common  practice  in  studying  fracture  of  the  unidirectional  composite 
to  assume  that  cracks  extend,  after  initiation,  along  a  direction  parallel  to 
the  major  material  axis,  i.e.  the  fiber  direction,  since  the  trans-fiber 
strength  of  the  composite  is  generally  much  lower  than  that  along  the  fiber 
direction.  Now  consider  the  problem  of  two  cracks  of  equal  length,  a/R  =  0.4, 
emanating  at  (R , 4>)  and  (R,ir+<f>)  from  the  edge  of  the  hole  in  the  glass /epoxy 
composite,  where  <j>  depends  upon  the  fiber  orientation  as  shown  in  Fig.  4. 
Fracture  parameters  such  as  mixed-mode  stress  intensity  factors,  K  and  K^, 
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which  describe  the  crack-tip  deformation  and  fracture  of  the  composite,  are 
calculated  explicitly  for  cases  of  various  fiber  orientations.  For  compara¬ 
tive  purposes,  corresponding  single-crack  problems  are  also  analyzed.  Results, 
expressed  in  terms  of  K.j./ (o ^/rraT)  and  K^/  ^Oco'^a ^  f°t  both  single  and  double 
crack  cases,  are  reported  in  Fig.  5.  The  associated  energy  release  rates  G 
are  determined  also  in  Fig.  6.  The  results  indicate  that  K^.  increases  contin¬ 
uously  with  the  increase  of  fiber  orientation  and  that  rises  to  a  maximum 

as  0  approaches  45°  and  then  gradually  decreases  to  zero.  The  opening  mode 
stress  intensity  factors  determined  for  most  cases  are  generally  larger  than 
that  of  the  shearing  mode.  For  the  particular  crack  length  studied,  differ¬ 
ences  in  Kj,  Kj.^-,  and  G  between  the  single  and  double  crack  cases  do  not  seem 
to  be  appreciable.  Nonetheless,  a  significant  discrepancy  may  be  expected  if 
the  a/R  ratio  becomes  larger,  which  will  be  discussed  in  the  next  section.  One 
extreme  situation  is  the  case  of  0  =  90  in  which  the  crack  estension  is  shown 
to  occur  along  the  direction  normal  to  the  loading,  i.e.  <f>  =  0°,  as  expected. 
This  gives  a  purely  opening  mode  fracture, i .e.  =  0.00,  due  to  the  symmetry 

of  loading  and  geometric  conditions. 

4.4  Crack  Extension  in  Off-Axis  Unidirectional  Composites 

Problems  concerning  crack  extension  from  the  edge  of  a  hole  in  the  com¬ 
posites  are  also  examined  in  detail  in  this  study.  The  case  of  an  off-axis 
unidirectional  glass/epoxy  with  fibers  oriented  along  9  —  45  is  presented 
here  for  an  illustrative  purpose.  The  geometry  and  loading  conditions  are  taken 
the  same  as  previous  examples  (Fig.  3) .  Both  single  and  double  crack  cases  are 
considered  again.  Solutions  for  the  problem  have  been  obtained  by  an  incre¬ 
mental  scheme  using  the  method  of  analysis  described  in  section  2.  Results 
given  in  Figs.  7  and  8  illustrate  the  fundamental  behavior  of  crack  extension 
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from  a  hole  in  the  off-axis  composites.  Effects  of  the  crack  length  on  stress 
intensity  solutions  and  energy  release  rates  are  shown  clearly  in  the  figures. 

Values  of  Kt/(/tT  a  )  remain  relatively  constant  for  a  wide  range  of  the  a/R 

I  00 

ratio  studied  and  change  significantly  for  both  smaller  and  larger  cracks  in 
the  single  crack  case.  For  the  double  crack  case,  K  /  (/rr"  oj  increases  mono- 
tonically  during  crack  extension.  It  is  interesting  to  note  that,  for  both 
single  and  double  crack  cases,  K  / (/ff  O  varies  approximately  linearly  with 
the  crack  length.  Slight  deviations  from  linearity  in  K  /(/m  cO  are  observed 
as  the  crack  (or  cracks)  extends  beyond  a/R  =1.  In  general,  the  results  indi¬ 
cate  that  all  K  ,  K^,  and  G  have  larger  values  in  a  double  crack  case  than  those 
in  the  corresponding  single  crack  problem,  presumably  due  to  the  larger  ’’effective1 
crack  length.  Moreover,  the  larger  the  crack  length,  the  greater  the  difference 
between  the  two  cases.  Composites  with  other  fiber  orientations  may 

exhibit  different  values  of  K's  and  G’s  during  crack  extension,  but  fundamental 
features  of  the  crack-tip  deformation  and  fracture  behavior  are  essentially 
the  same. 

Extension  of  this  method  to  more  complicated  fracture  problems  of  composites, 
such  as  viscoelastic  anisotropic  crack  analysis,  can  be  easily  achieved.  These 
further  developments  have  been  very  successful  and  are  reported  elsewhere  [24]. 
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5.  SUMMARY  AND  CONCLUSION 

An  analytical  investigation  of  cracks  emanating  from  a  circular  hole  in 
an  off-axis  unidirectional  fiber-reinforced  composite  has  been  conducted.  The 
method  of  analysis  employed  in  the  study  is  formulated  by  using  conservation 
laws  of  elasticity  and  fundamental  relationships  in  anisotropic  fracture 
mechanics.  The  current  approach  provides  a  very  convenient  and  accurate  means 
to  examine  the  complicated  crack  behavior,  when  used  in  conjunction  with  a 
suitable  numerical  scheme  such  as  the  finite  element  method.  The  problem  is 
eventually  reduced  to  a  system  of  linear  algebraic  equations  in  mixed-mode 
stress  intensity  factors.  Fracture  parameters,  describing  crack-tip  defor¬ 
mation  and  fracture  in  the  composite,  are  obtained  explicitly.  Effects  of 
material  anisotropy  and  crack/hole  geometry  are  examined  also.  Of  particular 
interest  are  the  energy  release  rates  associated  with  crack  extension;  their 
values  are  evlauated  for  various  cases. 

The  results  obtained  in  this  study  lead  to  the  following  conclusions: 

(1)  The  current  method  of  analysis  is  demonstrated  to  be  superior  to  other 
approaches  in  its  operational  simplicity  and  solution  accuracy  for  in¬ 
vestigating  mixed-mode  crack  problems  with  extremely  complex  crack  hole 
geometry  and  material  anisotropy. 

(2)  Mixed-mode  stress  intensity  factors  and  energy  release  rates  associated 
with  the  cracks  emanating  from  a  hole  change  very  appreciably  with  fiber 
orientation  in  the  composite.  K  and  G  increase  monotonically  with 
increasing  0;  but  K  reaches  its  maximum  at  0  =  45°,  and  then  decreases 
gradually  as  0  increases  further. 


17 


(3)  During  crack  extension  from  a  circular  hole  in  the  composite  with  a 
given  fiber  orientation,  the  opening  mode  stress  intensity  factors  are 
found  to  be  relatively  constant  for  a  certain  range  of  the  crack  length  (i.e. 

a/ R  <  1.0)  in  the  single  crack  case,  but  increase  rather  rapidly  in  the 
double  crack  case.  increases  linearly  with  the  crack  length,  but 

deviates  gradually  from  linearity  as  the  length  of  the  crack  exceeds  the 
hole  radius,  i.e  a/ R  >  1,  for  both  cases. 

(4)  The  energy  release  rates  and  mixed-mode  stress  intensity  solutions  for 
single  and  double  crack  cases  do  not  have  much  difference  for  problems 
with  small  crack-length-to-hole-radius  ratio  i.e.  a/R  <  0.5.  However, 
the  difference  increases  and  becomes  very  appreciable  as  a/R  increases 


during  crack  extension. 
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TABLE  1 


STRESS  INTENSITY  FACTORS  FOR  A  CRACK  EMANATING  FROM  A  CIRCULAR  HOLE 
IN  AN  ISOTROPIC  FINITE-DIMENSIONAL  PANEL* 


*4 

KII 

Stalk  and  Orringer  [20] 

2.302 

0.000 

Present  Solution 

2.333 

0.000 

*a/R  =  0 .4,  <t>  =  0° ,  56  eight-node  isoparametric  elements 
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